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We investigate the string breaking mechanism in «/ = 2 QCD. We discuss the lattice techniques 
used and present results on energy levels and mixing angle of the static BB\QQ two-state system. 
The string breaking is visualized, by means of an animation of the action density distribution as a 
function of the static colour source-antisource separation. 
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1. Introduction 

The breaking of the colour-electric string between two static sources is a prime example of a 
strong decay in QCD [j|]. Recently, we reported on an investigation of this two state system [Q, Bh, 
with a wave function \Q) created by a QQ operator and a wave function \B) created by a four-quark 
BB operator, where B = Qq. Q denotes a static source and q is a light quark. We determined the 
energy levels E\{r) — 2ra# and Eiir) — 2ms of the two physical eigenstates 1 1) and \2) which we 
decomposed into the components, 



|1) = cos0|0+sin0|B) 
\2) = -sin0|0+cos0|S). 



(1.1) 
(1.2) 



We characterize string breaking by the distance scale r c at which AE = Ej — E\ is minimized and by 
the energy gap AE C = AE(r c ). While these energy levels and r c are. first principles QCD predictions, 
the mixing angle 6 is (slightly) model dependent: within each (Fock) sector there are further radial 
and gluonic excitations and we truncated the basis after the four quark operator. 

We use rif = 2 Wilson fermions at a quark mass slightly smaller than the physical strange quark 
and a lattice spacing a = 0.166(2)ro ~ 0.083(1) fm and find [[J, r c = 2.5(3)ro 1.25(l)fm and 
AE C « 51(3) MeV, where the errors do not reflect the phenomenological uncertainty of assigning a 
physical scale to ro m 0.5 fm. Using data on the QQ potential and the static-light meson mass rag, 
obtained at different quark masses, we determine the real world estimate, r c = 1.13(10) (10) fm, 
where the errors reflect all systematics. An extrapolation of AE C however is impossible, without 
additional simulations at lighter quark masses. 

These results became possible by combining a variety of improvement techniques: the neces- 
sary all-to-all light quark propagators were calculated from the lowest eigenmodes of the Wilson- 
Dirac operator, multiplied by 75, after a variance reduced stochastic estimator correction step. The 
signal was improved by employing a fat link static action. Many off-axis distances were imple- 
mented to allow for a fine spatial resolution of the string breaking region. Last but not least great 
care was taken to optimize the overlap with the ground states, within the \ Q) and \B) sectors, using 
combinations of three dimensional APE and Wuppertal smearing. For details see Ref. We were 
able to achieve values of \a Q \ 2 = \(^ Q \Q)\ 2 = 0.62(2) and |a fi | 2 = K^elB)) 2 = 0.96(1) at r m r c 
for the overlaps of our test wave functions ^x) with the respective states on the right hand sides 
of Eqs. ( |l.l| ) and (L2). The almost optimal value of |ag| < 1 was essential to allow E\, £2 and 6 to 
be fitted from correlation matrix data, 
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(1.3) 



obtained at moderate Euclidean times: t >2a for Cbb and t > 4a for the remaining 2 matrix ele- 
ments Cqq and Cqb = Cbq at r w r c . 

In order to obtain dynamical information on the string breaking mechanism, we now study 
the spatial energy and action density distributions within the two state system. For instance one 
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can then ask questions about the localisation of the light qq pair that is created when r is increased 
beyond r c . The energy density will decrease fastest in those places where qq creation is most likely. 
Perturbation theory suggests that light pair creation close to one of the static sources is favoured by 
the Coulomb energy gain while aesthetic arguments might suggest a symmetric situation with qq 
dominantly being created near the centre. 

Another motivation for a more detailed study is our wish to relate the static limit results to 
strong decay rates of quarkonia. In the non-relativistic limit of heavy quarks, potential "models" 
provide us with the natural framework for such studies. In fact at short distances, r<l /A, potential 
"models" can be derived as an effective field theory, potential NRQCD, from QCD One can 
easily add a BB sector, as well as transition terms between the two sectors, to the QQ pNRQCD 
Lagrangian. Strong decays would then be a straightforward non-perturbative generalisation of 
the standard multipole treatment of radiative transitions in QED. Unfortunately, transitions such as 
T(45) — > BB can hardly be classed as "short distance" physics. So, some modelling is required. The 
natural starting point again is a two channel potential model which might still have some validity 
beyond the short distance regime. The transition rate would then be given by phase space times an 
overlap integral between the T(45) wave function and the BB continuum, where some additional 
input is required to constrain the sandwiched interaction term (for details see e.g. Ref. 

2. The Method 

We follow Ref. |Q] and define action and energy density distributions, 

<*(x) = +#„(«)], (2.1) 
^) = i[4(x)-f„(i)], (2.2) 

where 

* W = WA 2 (*)|„) - (A 2 ) = Jta - {A 2 } (23, 

We have suppressed the distance r from the above formulae and n = 1 denotes the ground state 
(dominantly QQ at r < r c ) and n = 2 the excitation (dominantly BB at r < r c ). Electric and magnetic 
fields are calculated from the plaquette, 

E 2 (x + lc£j = W 1 £ [u Xji4 + U x ^ m ] , (2.4) 

2/3 1 3 _ — _ — 

5 W = -TlE [VxiJ W + + V ^4,U + Vx-dt-4,ij\ , (2-5) 

a 4 ;=1 

where j = mod(/,3) + 1 and, 

/— — — t — t 

3 

We implemented two different operators with the same continuum limits: in one case we identified 
U x ,\l with the link U x ^ connecting x with x + afx. In addition we used smeared operators, 

= Psu(3) ( + £ U x , v U x+a ^ul afl v ) , (2.7) 

V J 



U x ,nv — -5-tr I U X ^U x+a fi^U x+a v ^U x+a j± +a y V j . (2.6) 
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where 7 = 0.4 and the sum is over all six staples, in the three forward and three backward directions. 
The 7- value was tuned to maximize the average plaquette, calculated from the smeared links. Pgu^) 
is a projection operator into the SU (3) group. For the un-smeared plaquette zo = 1 inEq. (2^6) while 
for smeared plaquettes zo = 1 + 0(a s ) is adjusted such that the vacuum expectation value of the 
average plaquette remains unchanged. 

The plaquette smearing enhances the signal/noise ratio. Due to this smearing and the fat link 
static action used, the peaks of the distributions around the source positions (that will diverge in 
the continuum limit) are less singular than in previous studies of SU (2) gauge theory at similar 
lattice spacings In the continuum limit the results from smeared and un-smeared plaquette 
probes will coincide, away from these self energy peaks. The draw back of plaquette smearing is 
that exact reflection positivity is violated. However, our wave functions are sufficiently optimized 
to compensate for this. 

We insert the E 2 (x,t) and B 2 (x,t) operators at position t/2 into the correlation matrix C(t), 
Eq. (13), see Ref. [||]. For evenlodd t /a- values we average E 2 \B 2 over the two adjacent time slices, 
respectively. Using the fitted ground state overlap ratio ag/ag and the mixing angle 8 as inputs, 
we calculate the action and energy density distributions Eqs. (2A) and (2.2) in the limit of large t 
via Eq. (13) from the measured matrix elements. The distributions agree within errors within the 
time range 3a < t < 6a. The results presented here are all based on our t = 4a analysis. 



3. Results 

To set the stage, we display the main results of Ref. ^ in Figure [Ij In the left figure we also 
speculate about the scenario in the real world with possible decays into BB as well as into B S B S . 
As discussed above, for our parameter settings and «/ = 2 string breaking occurs at a distance 
r c 1 .25 fm. In the right figure we show the mixing angle as a function of the distance. The BB 
content of the ground state is given by sin0. Within our statistical errors 6 reaches n/4 at r = r c . 
Remarkably, there is a significant four quark component in the ground state at r < r c while for 
r > r c the limit 6 — > n/2 is rapidly approached. 




r/fm 

Figure 1: The energy levels and the mixing angle in physical units for «/ = 2. The bands in the left figure 
reflect the expected nf = 2 + 1 scenario. 
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Figure 2: The transverse profile in the centre of the flux tube. 

The quality of our density distribution data is depicted in Figure ^| for the ground state at a 
separation slightly smaller than r c , as a function of the transverse distance x from the QQ axis. Due 
to cancellations between the magnetic and electric components the energy density is much smaller 
than the action density: for the comparison we have multiplied the energy density data by the 
arbitrary factor of 29. Note that the ratio o/e will diverge like — InaA in the continuum limit. The 
differences between the shapes of the energy and action density distributions are not statistically 
significant. Here we only visualize the more precise action density results. 

We employ several off-axis separations. Assuming rotational symmetry about the interquark 
axis, each point is labelled by two coordinates, x denotes the distance from the QQ axis and y 
denotes the longitudinal distance from the centre point. We define an interpolating rectangular grid 
with perpendicular lattice spacing a and the longitudinal spacing slightly scaled, such that the static 
sources always lie on integer grid coordinates. We then assign a quadratically interpolated value to 
each grid point z, obtained from points in the neighbourhood, |z — (x,y)\ < e = a. On the axis the 
data points are more sparse and we relax the condition to £ = \[%a while for the singular peaks we 
maintain the un-interpolated values. 

We append two mpeg animations of the action density distribution as a function of the distance 
lfm < r < 1.5 fm for the ground state (string fission) and the excited state (string fusion). The 
picture frames inbetween the measurement points of Figure [j] are linearly interpolated. Note that 
the time is not a linear function of the distance but dilated within the string breaking region. On 
a linear time scale string breaking takes place rather rapidly. The starting frames are displayed in 
Fig. |[ The colour encodes the relative statistical errors and the lattice mesh represents our spatial 
resolution a. In spite of the fact that the excited state has (within errors) the energy of two isolated 
static-light mesons, a significant string component is already present at r < r c which then grows as 
r is further increased. This is already obvious from the 6 (r) of Figure [j]. 
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Figure 3: Action density distribution for the ground state and the first excitation. 



4. Conclusion 

At any distance the action density looks like a superposition between (hypothetical) QQ and 
BB distributions with little interference: light pair creation seems to occur non-localized and in- 
stantaneously. Applying the decay model of Ref. where the interaction term is instantaneous 
and only depends on the separation r, we undershoot the experimental T(45') — > BB decay rate by a 
factor of about two. This appears very reasonable, considering the crudeness of the model and the 
fact that the gap AE C will increase with lighter, more realistic sea quark masses. We are studying 
the situation in more detail. 
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